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SUMMARY 


Pulsed Electrothermal (PET) propulsion is a technique for 
avoiding the material temperature limitations of steady 
electrothermal propu.'sion, by operating at megawatt power levels 
for microsecond time scales, with the time between pulses 
measured in milliseconds. By varying the energy per pulse and 
the pulse rate, operating parameters can always be found which 
prevent the electrodes from reaching the melting point. A second 
benefit of pulsed operation comes from the 1000 atmosphere 
pressure levels achieved during the pulse. The supersonic nozzle 
flow is then in an equilibrium rather than a frozen state, so 
that the ionization energy of the propellant is recovered in the 
nozzle. In addition, the high Reynolds number in the nozzle 
allows the nozzle flow to be treated as inviscid. 

This report is divided into two parts. Part one is an 
experimental study of an ablative wall PET thruster, using 
polyethylene propellant. In order to achieve high confidence in 
the performance measurements, a calibrator is used in conjunction 
with the thrust stand, and twin ti me-of-f 1 i ght probes are used to 
measure exhaust velocity. Energy delivered is measured from the 
time-integrated product of terminal current and voltage, and mass 
loss is determined with a highly accurate balance. The 
performance parameters are determinded by averaging the results 
for a ten shot series under constant conditions. 

The thruster is tested with four configurations of capillary 
length and pulse length. A 15 cm long capillary and 15 fi sec 
pulse are used to produce the unsteady mode condition, for which 
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the acoustic time in the capillary is greater than the pulse 
length tp (ctp/2i < 1). Quasi-steady conditions are created with 
4.5 cm and 5.5 cm capillaries, used with 48 //sec and 100 //sec 
pulses, for which ctp/2jt >> 1. 

As expected, the thruster runs at a relatively low average 
temperature in the unsteady mode, resulting in a 1000 sec 
specific impulse. The thrust to power ratio is .10 N/kW, and the 
thruster efficiency is .56. A large fraction of the exhaust mass 
in ths mode exists as "afterglow" in the tail of the pulse. 

Longer pul ses(quasi -steady mode) are tested to reduce the 
relative effect of the afterglow. The 48 //sec and 100 //sec 
pulses result in a 1500 to 1750 second specific impulse. The 
thrust to power ratio falls to .05 N/kW, and the efficiency 
averages .40. The exhaust velocity is measured as 27 km/sec, 
indicating that the afterglow effect is still important for the 
quasi-steady mode. 

Part II of the report addresses the propellant feed problem, 
and proposes a PET propulsion design for an o rb i t- ra i s i ng 
mission. A conceptual design of a 40 kW PET system is presented 
for a Space-Based Radar satellite. The selected mission is orbit 
raising of the 7000 kg SBR satellite in 15 days with a 
characteristic velocity of dV = 1000 m/sec. The point design 
system performance is .62 efficiency at 1000 seconds. The 
system's reliability is enhanced by incorporating 20 thruster 
modules of 20 kW each, fired in pairs. The thruster design is 
non-ablative, and uses water propellant from a central storage 
tank, injected through the cathode. System dry mass is 
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extrapolated from existing space-qualified hardware* and is 
estimated to be 10 kg/kw. 

I. INTRODUCTION 

In a previous report[l] we presented performance data for an 
ablative wall, pulsed el ectrotherma 1 (PET) thruster. In that 
device, a 2 kj electric discharge was used to create a 
polyethylene or Teflon plasma at a pressure of about 700 atm and 
a temperature of 2.5 eV. From current, voltage, impulse bit, 
exhaust velocity and mass loss measurements, we concluded that 
the PET thruster gives repeatable results, operates with 
equilibrium flow in the nozzle, and substantially recovers the 
plasma ionization energy. Based on calculations of losses from 
evaporation, dissociation, and unrecoverable thermal energy, 
these results demonstrated that the PET thruster has the 
potential of achieving efficiencies in the .70 to .80 range at 
2000 seconds specific impulse. 

The next steps in PET thruster development are discussed in 
this report. In the first reportfl], it was concluded that 
refinements in the performance measurement techniques could be 
easily made that would increase accuracy and confidence in the 
results. We have therefore developed a precise calibration 
system that simulates the impulse bit with an overall accuracy of 
better than 1 %. We have also used twin photodiode probes to 
improve the exhaust velocity measurement, and have improved the 
thrust stand design to reduce mechanical noise on the thrust 
signal. Improvements have also been made in the technique for 
measuring mass loss. 
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The second important step for PET thruster development, 
addressed in this report, is to solve the propellant feed 
problem. Although solid propellant feed systems for the Teflon 
Pulsed Plasma Thruster have been successfully developed and 
orbited, it was not at all clear that this could be easily done 
for the PET thruster's coaxial geometry. As a focus for a feed 
system conceptual design effort, it was decided to select a 
representative mission, for which the PET thruster is 
particularly attractive. The result of this study is a 40 kW PET 
propulsion system for the Lincoln Labs' Space Based Radar, using 
water propellant. The results of the experimental effort and the 
conceptual design are presented below and in Appendix A. 


II. PET THRUSTER PERFORMANCE MEASUREMENTS 


In a previous report[l], diagnostic techniques for PET 
thruster performance measurements were discussed in detail. This 
discussion is summarized in Appendix A of this report. We have 
since developed several techniques for refining the PET thruster 
performance measurements. 

The performance of the PET thruster is defined through five 
primary diagnostic techniques. These are: 

1. Discharge current, measured by an RC integrated Rogowski 
loop. 

2. Discharge voltage, measured by a resistive 935:1 
voltage divider. 

3. Impulse bit, measured by the recoil momentum of the 
thrust stand. 

4. Exhaust velocity, measured by ti me-of-f 1 ight to three 
types of probes at different positions in the exhaust 
stream. 

5. Mass loss, determined by accurate weighing of the 
discharge liner, cathode, and nozzle. 

Typical profiles for current, voltage. Impulse bit, and 
exhaust velocity probes are shown in Figures 1-3. The analysis 
of PET thruster performance is based on terminal measurements. 
The current and voltage are multiplied to give input power (Fig. 
1) and then integrated to give the delivered energy, Ep. 
Typically, the discharge resistance is at least ten times larger 
than the external resistance, so that E D / ( 1 / 2 C V 0 2 ) > 0.9. The 
slope of the recoil signal (Fig. 2) is determined graphically, 
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PET THRUSTER TERMINAL MEASUREMENTS 
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and is multiplied by the recoil mass to give impulse bit. These 
two measurements immediately give the thrust to power ratio T/P 


from 


V WW * W 


Previously reported results for polyethylene propellant have 
given values of T/P ranging from .0 5 [N/kW] to .10 [N/kW] 
depending on the pulse length. 

The exhaust velocity and mass loss measurements are required 
to determine specific impulse and thruster efficiency. Their 
product should be equal to the measured impulse bit: 

/Tdt = Amu g 

As was shown previously, however, the quantity 4mu e exceeds the 
impulse bit /Tdt by nearly a factor of 2. The explanation for 
this discrepancy is that the mass velocity measured by the time- 
of-flight exhaust probes is higher than the mean mass velocity 
for the entire pulse. In support of this, a computational model 
of the discharge shows that about two-thirds of the mass is 
ejected at high temperature during the pulse, and one third is 
ejected at much lower temperature and velocity in the pulse 
"afterglow" (Fig. 4 ). 

Since the afterglow portion of the pulse reduces the mean 
exhaust velocity, it was decided to minimize the importance of 
the afterglow by increasing pulse length. We report here on two 
test series, of ten shots each, using polyethylene propellar-t. 
One series was run with a 48 fisec pulse, duplicating the 48 n sec 
tests reported earl1er[l]. The second series was run with the 
same capillary geometry but with a 100 /isec pulse. The results 
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COMPUTED CAPILLARY MASS LOSS 



Figure 








of these tests are summarized below, and are reported in more 
detail in Appendix A. 

Since a major goal of the present experiments is to improve 
measurement accuracy, new approaches were taken to the techniques 
for measuring thrust, exhaust velocity and ablated mass. These 
are discussed below. 

Thrust 

Since the thrust measurement is a vital performance 
parameter, a means of calibrating the thrust st*nd has been 
devised (Figs. 5 and 6). A physical pendulum of mass m is 
constructed to rotate about pivot point 0, and impact the 
thruster, mass M, at a radius r. The impact point at r is 
designed to be the center of percussion about 0 so that no impact 
forces are transferred to the suspension structure. The pendulum 
and thrust stand are then considered a purely rotational system, 
with angular momentum about 0 being conserved. 

The pendulum velocity immediately before and after impact is 
found by passing the focussed beam of a helium-neon laser through 
a grid of equ* -spaced slits (50 line/ inch Edmund Scientific 
Ronchi Ruling), and measuring the chopping frequency of the 
transmitted beam (Fig. 6). All other fixed character!' sti cs of 
the pendulum such as the moment of inertia are known, so that the 
net momentum transferred to the thrust stand can be calculated. 
Simultaneous - ^, the recoil momentum of the thrust stand is found 
using the position transducer (Fig. 7). Comparison of the two 
measurements yields a calibration constant that is applied to all 
subsequent thrust data obtained using the transducer method. 





CALIBRATOR OPTICAL SCHEMATIC 



Figure 6 



The design goal for the calibrator was to establish absolute 
confidence in the thrust measurement at the level of 5% or 
better. Some features of the calibration system are therefore 
noteworthy. All the basic physical characteristics of the 
pendulum, such as the mass m, moment of inertia I, and the 
various distances indicated in Figure 5, were measured to about 

0. 1% accuracy. Table I summarizes the relevant data. In Table 

1, T is the period of the free swinging pendulum, and Ax 
characterizes the Ronchi ruling grid pattern. The incident and 
rebound swing velocities, v^ and Vg , are determined using the 
known Ax and by a linear least squares fit to the timing data 
presented by the rising and falling edges of the chopped laser 
response (Fig. 7). Overall, the pendulum calibrator has an 
inherent accuracy of ± 2%. 

In order to mimic the recoil experienced by the thrust stand 
during a discharge, the pendulum was arranged to strike directly 
on and parallel to the thruster axis. As well, both the pendulum 
and thruster have hardened steel inserts at the contact point. 
Independent measurements indicate that this results in an impact 
time that is comparable to the duration of the discharge current 
pul se . 

Throughout the course of the thruster experiments, 
calibration runs were made in order to monitor the stability and 
accuracy of the thrust measurements. The calculated recoil 
velocity differed from that measured by the position transducer 
by ±3% of the mean in the worst case. The average difference 
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Figure 7 




Table I. Pendulum Characteristics 


g = 980.1 cm s“^ 

± 0.10% 

h = 12.09 cm 

± 0.41% 

m * 404.81 <jm 

± 0.05% 

T = 0.7875 s 

± 0.08% 

M = 2840 gm 

± 0.18% 

r = 15.24 cm 

± 0.28% 

H = 19.77 cm 

± 0.13% 

Ax = 254.0 yin 

± 0.8% 


V r = MrH ^ V 1 
I * n,gh 
V-J = AX/At.j 
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over several calibrations was only ±1.3%. 

A further improvement in measuring the thruster recoil 
velocity was made by simply relocating the position transducer 
and target plate. Previously, experiments were performed by 
having the position transducer 'view 1 the recoil motion directly 
along the thrust axis (see Fig. 4 of Appendix A). In such a 
configuration, axial acoustic modes of the thrust stand produced 
a large amplitude, high frequency oscillation superimposed on the 
relatively slowly varying recoil signal. (See, for example. 
Figures 7 and 8 of Ref. [1].) However, by locating the 
transducer and target as shown in Figure 5, this acoustic 'noise' 
has been significantly reduced, allowing the recoil momentum to 
be determined to better than ±4% (Fig. 2). 

Exhaust Velocity 

In order to measure nozzle exhaust velocity, three time-of- 
flight techniques were developed: ion saturation probes, piezo- 
electric stagnation pressure probes, and photodiode probes. The 
exhaust velocities measured using each of these techniques are in 
good agreement, and a detailed description of the individual 
methods is given in Appendix A. 

The primary velocity diagnostic selected was the photodiode 
probe, since this type of probe produced the most reliable and 
reproducible time-of-arri val signals. Two identical probes were 
installed downstream of the nozzle and separated axially by 50.0 
cm. By using two probes simultaneously and noting the difference 
in the arrival time of the exhaust plasma, the exhaust velocity 
is determined for each shot in the test series. As an example. 
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the time difference At and exhaust velocity for each shot in the 
long pulse (100 //sec) test series is shown in T able II. The 
major uncertainty here lies in extracting At from the individual 
photodiode traces. Each value for At is accurate to < ± 10% so 
that, over the 10 shot test series, the average exhaust velocity 
is established to ± 5%. The tank pressure was kept at less than 
0.003 torr in all cases, to prevent interference by the 
background gas with the exhaust gas. 

Mass Loss Measurements 

Ablation studies are carried out via a simple weighing 
procedure. All components exposed directly to the discharge 
plasma (cathode, capillary and nozzle) are weighed before and 
after each series of test shots. The improvement made here is to 
include the nozzle in the mass loss data, since the nozzle 
structure incorporates a polyethylene liner in the vicinity of 
the throat, where some ablation might still be expected. 

All mass measurements are made using a Voland Model 220R 
balance having ± 0.1 mg resolution. However, experience has 
shown that the thruster assembly, disassembly and cleaning 
procedures introduce a total uncertainty of 1-2 mg. Therefore, 
capillary loss is typically accurate to < ± 5% whereas the 

cathode and nozzle erosion data is less precise due to the 
smaller total mass loss. 

The performance of the PET thruster was evaluated in 
four series of ten shots each. The capillary length and current 
pulse length were varied to permit the thruster to run in an 
unsteady or quasi-steady mode, depending on the number of 
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Table II. Exhaust Velocity Data For 100 ys 
Pulse Test Series 


Shot No. 

u £ ( km/sec) 

At(yS) 

P165 

- 

- 

166 

27.8 

18.0 

167 

27.9 

17.9 

168 

27.6 

18.1 

169 

24.8 

20.2 

170 

25.4 

19.7 

171 

26.2 

19.1 

172 

25.6 

19.5 

173 

28.1 

17.8 

174 

29.1 

17.2 


< u e > =26.9 km/sec 
a = ± 1.47 km/sec ± 5.5% 


Table III. Experimental Parameters for 
Polyethylene Propellant 


Test 

Series 

Capillary 
Energy, kJ 

Pulse 

Length, ys 

Capillary 
Length, cm 

ct 

1 . 

1.80 

15 

15 

0.4 

Unsteady 

2. 

1.78 

48 

5.5 

3.5 

Quasi- 

steady 

3. 

1.89 

48 

'•.5 

4.3 

Quasi- 

steady 

4. 

3.98 

100 

4.5 

8.9 


Quasi- 

steady 
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acoustic travel times ct p /2j 2 during the current pulse. The four 
test series are summarized in Table III below. 

After each series, the results of the ten shots were 
averaged. The performance results are summarized below in Table 
IV. 

Based on the experimental results, the estimated precision 
of the results in Table IV is: 

1. Recoil motion x(t); impulse bit is accurate to 5%. 

2. Exhaust probe time-of-f 1 i ght ; exhaust velocity is 
accurate to 10%. 

3. Current; accurate to 3%. 

4. Voltage; accurate to 3%. 

5. Mass loss 4m; accurate to 6%. 

Summary 

The probe measurements of exhaust velocity give results 
which are consistent with our understanding of the physics of the 
discharge. In the unsteady mode case (Series 1), the current 
impulsively heats the discharge plasma, which is then accelerated 
as the temperature decays. The measured exhaust velocities are 
20-21 km/sec. The quasi-steady cases (Series 2-4) achieve nearly 
steady state conditions, and hence a higher mean sound speed for 
the duration of the pulse. The exhaust velocities are 
correspondingly higher: 26-28 km/sec. A discussion of the 

agreement of these measured velocities with a computational model 
is given in Appendix A. 

The thrust to power ratio [N/kW] can be evaluated in the 

form 
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Table IV. PET Thruster Performance With Polyethylene, 10 Shot Averages 
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The values of T/P for the test shots (Table IV) are shown below 

In Table V. The velocities shown are these measured by the 

velocity probes for the mass exhausted during the current pulse. 

The thruster efficiency r» j 1 s conventionally defined as 

_ exhaust kinetic energy _ mu e 
n T energy delivered ~ Ep 

Table V. Measured Thrust to Power Ratio 


for Polyethylene 




u 0 ( Probe) 

T/P 

1 . 

Unsteady, 

15 psec pulse 

21 km/sec 

.10 N/kW 

2. 

Quasi-steady 
48 usee pulse 

27 km/ sec 

.053 N/kW 

3. 

Quasi -steady 
48 psec pulse 

26 km/ sec 

.046 N/kW 

4. 

Quasi -steady 
100 psec pulse 

27 km/ sec 

.048 N/kW 


For a pulsed device, the exhaust kinetic energy is 

1/2 mu | = [imu|]dt 


where the integral is continued past the end of the current pulse 
to account for the discharge afterglow, during which propellant 
evaporation decays to zero. Unfortunately the time dependence of 
m and u are not known, so that we must define the efficiency in 
terms of impulse bit as: 


- C/Tdt] 
'T " 2amE 
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w‘;c re 4m is the ca p i 1 1 a ry ablated * ' s , and is used to define the 

spe c impel se by ^ 

*sp gAn~ 

The :=sta from lable V above c<r cerefore be used to evaluate 
n a r, d I for Series 1-*, >.iown below in Table VI. It is 

important to note that we basing the efficiency and I S p in 
Table VI on the ablated sna.s and not on the measured exhaust 
velocity. 

The specific impulse values in Table VI are about half of 
those indicated by the exhaust probe velocity measurements. This 
is primarily because the mass fraction ejected in the afterglow 
tail of the pulse has significantly lower velocity than the 
fraction ejected during the pulse. This is expected since the 

afterglow plasma vs cooler. 


Table VI. PET Thruster Efficiency and 
I sp for Polyethylene at 2 kJ 




n T 

*sp 

1 . 

Unsteady 
15 usee pulse 

.56 

1000 sec 

2. 

Quasi -s tea y 
48 usee pulse 

.39 

1500 sec 

3. 

Quasi -steady 
48 usee pulse 

.37 

1600 sec 

4. 

Ojas'* -steady 
100 usee pulse 

.42 

1750 sec 
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Concl usi ons on PET Thruster Perf ormarce Measurements 
A single-shot PET thruster has been built and tested on a 
thrust stand, giving repeatable reuits over a series of shots. 
Exhaust velocities for polyethylene propellant vary from 20 to 27 
hii'/sec, as measured by exhaust probes. These velocities are 
consistent with computational predictions of capillary conditions 
with an adiabatic expansion to high Mach number in a conical 
nozzle. Single pulse spcific impulse nd efficiency 
measurements, based on ablated mass, show a thruster efficiency 
of .37 to .56 in the 1000-1750 second range- We believe that an 
understanding of loss mechanisms in the discharge would lead to 
an improved design with a thruster efficiency of .70 to .80. 


III. PET THRUSTER SYSTEM CONCEPTUAL DESIGN 

In order to evaluate the PET thruster concept as a system, 
we have performed a conceptual design for a PET propulsion 
mission [2]. The two principal goals of the design are one, to 
design a practical propellant feed system and two, to estimate 
the system component masses. 

The selected mission for the conceptual design is the orbit 
raising of a Space-Based Radar Satellite, a concept being 
developed by MIT Lincoln Laboratories (Fig. 8) [3]. The mission 
for this design is to raise the orbit of the SBR and change the 
crbital plane, following deployment by the Space Shuttle. The 
orbital altitude gain is from 150 n.m. circular to 600 n.m. 
circular, and the plane change is 4°. For the purposes of this 
study, the mission characteri Stic velocity is taken as 4V = 1000 
m/sec . 

Several aspects of this SBR orbit-raising mission make it 
appropriate for the PET thruster conceptual design. 

1. Solar panels and batteries on the SBR (Fig. 8) can 
supply 50 kW of power on a 200V DC bus. This power can be made 
available for orbit raising propulsion. 

2. As much as 30 days can be allocated for the orbit 
raising, well within the capability of PET propulsion. 

3. The proposed SBR initial orbit (150 n.m. circular, 57 
inclination) is difficult to reach by the Shuttle, placing a 
premium on initial mass. The low propellant mass of electric 
propulsion is therefore attractive. 







4. The inherently low acceleration of a PET-propel led SBR 
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system, approxi mately 10 g, is consistent with the structural 
requirements of a large deployed solar array. 

5. The high specific impulse and pulsed operational 
capability of PET propulsion would be attractive for station- 
keeping and precision attitude control during radar operations, 
although this aspect was not included in this design. 

6. PET propellants can be benign chemically, an important 
consideration with respect to contamination of solar panels by 
propulsion system exhaust. 

Mission Parameters 

A Lincoln Laboratory estimate of the mass of the SBR 
satellite with chemical propulsion is 9800 kg, of which 2800 kg 
is propellant. These values are consistent with a characteri Stic 
velocity of AV = 1000 m/sec and a specific impulse of 300 
seconds. By comparison, the specific impulse of PET is 
considerably higher, resulting in increased trip time and reduced 
propellant mass for the orbit-raising mission. 

The calculation of trip time for the SBR satellite is 
performed for the simple case of constant specific impulse, 
thrust, and efficiency. The SBR satellite has sufficient battery 
capacity so that thrust can be maintained during solar eclipse. 
Hence, trip time and thrust time are assumed equal. The thrust 
is given by 2n c n T P 


where n c and n-y are the power conditioning and PET thruster 



efficiencies, P is available bus power, and u g is exhaust 
velocity. For the simple case of constant T and u 0 , the 
integration of T = m(t)a gives a trip time of At of: 


At[sec] = 


1/2 m f u| 




- 1 ] 


The trip time and required initial mass m.. are calculated for the 
values given below in Table VII. The calculation of the PET 
propulsion system dry mass of 400 kg is discussed in detail 
bel ow. 


Table VII. Space Based Radar Trip 
Time Parameters 

Final mass, m f 7400 kg 

(SBR, 7000 kg; PET 
propulsion, 400 kg) 

Delivered Power, 40 kW 


Thruster efficiency, n T 0.5, 0.6, 0.7 


Mission aV 
Specific impulse 


1 km/sec 

400 <1 <1600 sec 
sp 
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The trip time for the &i = 1000 m/sec orbit transfer is 
plotted in Fig. 9 for the parameters in Table VII. The thruster 
efficiency range of .50 to .70 represents the expected operating 
regime of the PET thruster. Trip time in Fig. 9 decreases 
mo n o to n i c a 1 1 y with decreasing I sp , reaching 7-9 days for 
continuous thrust at 400 seconds, for an initial mass of 9500 kq. 
At higher specific impulse tne initial mass is reduced, reaching 
8250 kg at 1000 seconds. At still higher I sp the initial mass 
decreases further but the trip time increases rapidly. A design 
point is selected giving At = 15 days at 1000 seconds and .70 
efficiency (Fig. 9). Recent experimental results, discussed 
above, indicate that this design point ;s a realistic one for a 
PET propulsion system. A summary of the design point conditions 
is shown below in Table VIII. 

Propul sion System Desi gn 

Experimental work to date on the Pulsed Electrothermal (PET) 
thruster has been limited to ablative-wall versions of the 
device, using polyethylene or Teflon propellant, and described 
previously [1,4,5]. For the SBR orbit-raising mission the 

thruster must be non-ablative to meet the required system life of 
15 days. Furthermore, a liquid propellant is required so that a 
central tank system can be used to supply several thrusters. 
Since the thrusters are low mass and compact, the use of liquid 
propellant permits several tens of thrusters to be installed, 
although only two are operated at any given time. This approach 
greatly increases system reliability with a low mass penalty. 
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Table VIII. Design Point Parameters for 
SBR Mission, aV = 1 km/sec 


Initial mass, m.. 

8250 kg 

Final mass, m t 

7400 kg 

Propellant mass, m p 

850 kg 

Thruster efficiency, n T 

0.70 

System efficiency. 

0.62 

n = n c n T 


Specific impulse, I 

1000 sec 

Input power, P 

45 kW 

Thrust, T 

5.6 N 



Water (mol. wt. 18) Is the liquid of choice for the 
propellant. The advantages of water are: 

1. Easily storable. 

2. Existing space-qualified flow hardware. 

3. Simple blowdown injection. 

4. Propellant chemistry matched to mission. 

5. Low expected spacecraf t/envi ronmental contamination. 

6. Non-reactive, non-toxic. 

7. Convenient vapor pressure for no-switch operation. 

8. Known thermodynamic properties. 

Items 4 and 7 in the above list require further comment. 
First (item 4), the propellant chemistry of water is well matched 
to the SBR mission, which is designed for 1000 seconds specific 
impulse. This specific impulse can be achieved by water at a 
temperature below 1 eV (11,6000 K), as required for non-ablative 
thruster operation. Unsteady heat transfer calculations for ten 
microsecond, 1 eV discharges show that the working temperature 
will not be exceeded for available insulator, electrode and 
nozzle materials. 




Numerical calculations of the capillary performance with 
water have been performed (Fig. 10), and water can be seen to be 
similar to the behavior of polyethylene in terms of pressure, 
temperature and sound speed. All curves in Figure 10 were 
calculated using the experimental power pulse shown at the bottom 
of Fi gu re 1 . 

Second (item 7), the PET thruster is designed to fire at a 
rate of several hundred pulses per second. This is accomplished 
without the aid of a switch or trigger circuit, by permitting 
Paschen breakdown through the water vapor surrounding the 
injected liquid. This mode of operation has been experimentally 
demonstrated recently at GT-Devices at breakdown voltages of 
several kilovolts. 

The basic elements of the PET thruster operating with water 
are shown in Figures 11-14. In the schematic of the operating 
system (Fig. 11), a high voltage DC charging supply, powered from 
the SBR satellite 200V DC bus, charges a capacitive pulse forming 
network. Simultaneously, liquid water flows into the discharge 
chamber through the central electrode. The flow rate and the 
charge rate are roughly timed, so that the breakdown voltage is 
reached when the leading edge of the injected water is just 
reaching the throat of the supersonic nozzle. The PFN energy is 
discharged in a few tens of microseconds, and the heated water 
expands through the equilibrium flow nozzle. The cycle then 
repeats. Because no pulsed water valve is in the system, the 
pulse rate can be several hundred pulses per second. 
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Figure 11 






The water-injected PET thruster requires a non-ablative 
Insulator capable of withstanding thermal and mechanical shock, 
and the high heat transfer conditions of the discharge. The best 
candidates for this insulator are high-performance ceramics. 
Hot-pressed silicon nitride, for example, has been used for gas 
turbine wheels and diesel pistons [6], and is one of several good 
candidates for the PET thruster application. Advantage can be 
taken of the very high compressive strength of ceramics by 
shrinking a steel compression jacket (Fig. 11) around the 
insulator, as has been done in similar applications [7]. 

A schematic of the PET thruster module with a radiator is 
shown in Fig. 12. Two modules are operated in parallel in the 
present design. A high voltage insulator is used to isolate the 
water feed system from the high voltage current pulse. This 
insulator takes advantage of the high dielectric strength and 
resistivity of water. 

The mass of the SBR 40 kW PET propulsion system is estimated 
from the component masses of space qualified hardware. We have 
adapted a design similar to that for the Viking 75 Orbiter 
spacecraft, which carried a helium-pressurized chemical 
bipropellant system for trajectory correction and maneuvering in 
orbit around Mars [8]. Total Viking Orbiter propellant mass is 
1426 kg, and the mass of propellant and pressurant tanks for the 
850 kg propellant PET system is scaled down proportl onat^ly . 
This mass estimate is conservative, since the V0'75 design 
incorporates two tanks with an internal pressure of 8 atm, 
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whereas the PET system requires one tank with an interna] 
pressure of 4 atm. 

The SBR 40 kW PET design achieves redundancy by 
i ncorporati ng 20 PET thrusters with propellant shutc f f valves and 
electrical connect switches. Two 20 kW thrusters are in 
operation at any given time, each driven by an independent 100 J 
pulse forming network (Fig. 13). 

The flow and electrical schematics are shown in Fig. 14. 
The pressurant helium for the propellt tank is regulated to 4 
atm, and the pressurized water to each thruster flows through a 
filter, a shutoff valve, and the injection orifice to the 
discharge chamber. Back flow caused by chamber discharge 
pressure is not expected to be significant, due to the small size 
of the injection orifice (<0.3 mm dia), short duration o f the 
pulse (10 //sec), and low density of the mass in the discharge 

3 

chamber (.005 g/cm ). Mass estimates for propellant flow 
hardware are based on the V0'75 system (19 items, total mass 12 
kg) [8] and on mass values supplied by manufacturers for 

solenoid valves and filters [9], cartridge-actuated valves [10], 
and high voltage switches [11]. 

Each thruster is connected to one of two 100 J pulse forming 
networks through a latching connect switch. Switch operation is 
only permitted with the PFN at zero voltage to prevent electrode 
arcing. In addition, each capacitor and inductor in the PFN are 
connected via DPDT latching switches, to allow removal of a 
failed capacitor from the system. Capacitor mass is based on a 
conservative 12 J/kg, which should permit a design life of 10^ 
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cycles in a non-reversing application [12]. This value can be 
compared with the capacitor bank used for the Pulsed Plasma 
leflon thruster which operates at 64 J/kg [13]. This 
conservative 12 J/kg rating is feasible because of the low energy 
per pulse of the PET system. 

The mass of the power conditioning and logic circuits is an 
extrapolation of series resonant charging technology. The SBR DC 
bus voltage is 200V, and the discharge breakdown voltage is 4 kV. 
It is estimated, from experience with Pulsed Plasma thruster 
technology, that an 89% efficient, 20 kW power conditioning and 
logic module can be built at a specific mass of 2.5 kg/kW [14]. 
The power conditioner is thus the most massive component in the 
system. 

3 

The thruster modules are designed tor containment of the 10 

O 

atm discharge pressure in a 1 cm capillary volume. They are 
conservatively estimated to have a mass of 1 kg each, including 
nozzles and insulators. The thrusters are mounted on a common 
radiator ring designed to radiate a total of 4 kW from 2 
thrusters at a mean temperature of 540°C. The radiator material 
is assumed to be a copper- s ta i n 1 e s s steel sandwich for high 
temperature strength, formed into a dished donut (Fig. 13). 

The mass of the structure is estimated from the V0'75 
design, for which structure mass was 11% of the total non- 
structural non-fluid mass [8]. A figure of 11% is therefore 
adopted for this design. 

The component masses for the SBR 40 kW PET propulsion system 
are summarized below in Table IX. The predicted specific mass of 




Table IX. Component Masses, SBR 40 kW 
PET Propulsion 


Mass, kg 

A. Propellant Handling 


Water tank, 860 kg, 4 atm 

61 

Water holdup 

10 

Helium pressurant tank, 230 atm 

21 

Helium pressurant 

3 

Main shutoff, regulator, relief, 

6 

main filter 


Thruster filters, 20 

10 

Thruster shutoff valves, 20 

10 

Propel 1 arc lines 

10 

Subtotal 

131 kg 

Power Conditioning 


Series resonant DC-DC charger, 

100 

2 § 20 kW, 5 kV 


Capacitors, 20 @ 10 J 

16 

Inductors, wiring, insulation 

10 

Thruster connect latching relays, 20 

30 

PFN connect switches latching DPPT, 20 

25 

Subtotal 

181 kg 

Thrusters 


Modules with nozzles, 20 

20 

Radiator ring, 20 thrusters 

25 

Subtotal 

45 kg 

Structure, 11% of non-fluid mass 

40 kg 

Total mass: 

400 kg 

Specific mass: a 

= 10 kg/kW 
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the entire system is a = 10 kg/kW. 

If more thruster redundancy Is required for longer syste „ 
life, 20 kW thruster modules can be added at a total system mass 
penalty of 7 kG per thruster. Doubling the number of thrusters 
from 20 to 40, for example, would increase the total mass from 
400 kg to 540 kg and raise the system specific mass to 13.5 
kg/kW. 

Of the major research issues which need to be addressed for 
PET propulsion thruster efficiency is the only one which can 
have a large impact on the design. The remaining issues, such as 
peak discharge temperature and pressure, water injection nozzle 
design, and insulator and capacitor life, will have a relatively 
small effect. If the present experimental efficiency of .50 
cannot be improved to the assumed value of .70, the trip time 
would increase from 15 days to 24 days, but there would be no 
increase in the initial mass. 
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IV. CONCLUSIONS 



We have performed definitive and accurate experiments to 
characterize the performance of a single-shot, ablative Pulsed 
Electrothermal (PET) thruster. The data show that the thruster 
efficiency range is .37 to .56 at 1000 to 1750 seconds. In our 
opinion, the efficiency can be raised to .70 through design 
improvements . 

We have also proposed a non-ablative PET thruster concept, 
using water propellant. Numerical calculations show that this 
thruster would have similar performance to the ablative thruster, 
with significant system advantages in redundancy, reliability, 
and thruster lifetime. Based on these results, a repetitively 
operated, non-ablacive PET thruster can now be designed and 
tested, with a performance goal of .70 efficiency at 1000 
seconds. We expect that this device would prove superior to 
other propulsion devices operating at this specific impulse. 
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Ab stract 

Exhaust velocity, thrust, and mass 
consumption measurements are performed on a 
pulsed Electrothermal (PET) thruster using 
polyethylene propellant. The data verify 
theoretical predictions of equilibrium flow 
in the supersonic nozzle, resulting in 
nearly complete recovery of the energy of 
ionization. The thruster is tested in an 
acoustically unsteady mode (15 ysec current 
pulse) and in an acoustically quasi-steady 
mode (48 and 100 psec current pulses) at 2 
and A kj. The exhaust velocity of the 
propellant mass exiting during the current 
pulse is measured with three types of tlme- 
of-flight probes, and the impulse bit Is 
measured on a thrust stand. The measured 
thrust to power ratio is T/P - 0.10 N/kW at 
21 km/sec in the unsteady mode and T/P • 
0.05 N/kW at 27 km/sec in the quasi-steady 
mode. The discharge pressure, temperature, 
and sound speed predicted by a computation- 
al model are consistent with the measured 
thrust and discharge resistance. Propel- 
lant consumption is determined by weighing 
the loss for 10 shots at fixed conditions. 
The total measured consumption is in excel- 
lent agreement with that given by the mod- 
el, which predicts that about 2/3 of the 
mass is ablated during the current pulse 
and about 1/3 at later times. The computed 
sound speed for polyethylene is 9 km/sec, 
indicating that the exhaust plasma (27 
km/sec) is highly supersonic. The thruster 
efficiency varies from .37 to .56 for 1000- 
1750 ssconds specific impulse. 

Introduction 

In a previous paper [1] we have dis- 
cussed the theory of the Pulsed Electro- 
thermal (PET) thruster, a device which in 
principle can operate with 70* efficiency 
at a specific Impulse of 1000 seconds and 
higher. This level of performance would be 
particularly attractive for orbit raising 
of large satellites and other near-earth 
missions, which cannot be easily accom- 
plished by chemical propulsion. We have 
reoently completed a series of experiments 
to measure the single-shot performance of a 
PET thruster on a thrust stand. These 
experiments confirm that our theoretical 
understanding of the pulsed electrothermal 
approach is essentially correct, and that 
the PET thruster is potentially capable of 
significant performance improvements over 
steady-state devices. 
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In its simplest form, the PET thruster 
(Fig. 1) consists of a long, thin, dis- 
charge channel, closed on one end and cou- 
pled on the other end to* a conventional 
supersonic nozzle. The current flows from 
the a/iode located in the nozzle downstream 
of the throat, to the cathode located in 
the closed end. The pressure chamber inner 
wall is lined with an easily ablated, low 
molecular weight material such as polyeth- 
ylene or Teflon. The pressure chamber is a 
high strength dielectric, capable of with- 
standing plasma pressure pulses up to sev- 
eral thousand atmospheres. The high pres- 
sure develops a jet of high velocity plasma 
from the nozzle, which provides rocket 
thrust . 

The electrical energy store (Fig. 1) 
consists of a high voltage capacitive pulse 
forming network, capable of storing a few 
kilojoules at several kilovolts. The im- 
pedance of this source is matched to the 
discharge resistance R for efficient coupl- 
ing. The switch is an ignitron with resis- 
tance R <<R , so that over 90* of the stored 
8 

energy is transferred to the discharge in a 
single pulse. The width of the current 
pulse is 15-100 ysec, with a peak in the 
15-22kA range. Terminal voltage Is 2-3 kV, 
giving a peak power dissipation on the 
order of 50 MW. 

We discuss results from two PET thrust- 
er operating modes in this paper. The mode 
is determined by comparing the length of 
the current pulse, t , to the two-way lon- 
gitudinal acoustic Eime in the capillary, 
2 i/c . In the short pulse or unsteady mode, 

t c/2£<1, and steady flow is not achieved 
through the nozzle. In the long pulse or 
quasi-steady mode, t c/21»1, and steady 
flow is achieved during the high current 
portion of the pulse. The three current 
waveforms tested, with pulse lengths of 15, 
48 and 100 nsec, are shown in Fig. 2. 


Physios of Operation 


For a PET thruster operating with an 
ablative liner at millltorr background 
pressures, the initial voltage breakdown 
occurs as a surface flashover, mostly In 
carbon residue along the wall. Within 
several mloroseoonds the current, power 
input, pressure and temperature all rise to 
their peak values, and wall material begins 
to boil off, filling the oaplllary. The 
pressure, on the order of 10 5 atmospheres, 
generates a mass flow through the nozzle, 
which for constant power input balances the 
mass evaporating off the wall. When the 
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current falls at the end of the pulse, 
power input falls rapidly. The remaining 
plasma energy causes an afterglow of con- 
tinued wall ablation at low pressure for 
many microseconds, and this afterglow in- 
creases the total mass ablated during the 
pulse . 

Because the currents in the discharge 
are typically 20kA, the magnetic pressure 
is only a few percent of the gasdynamic 

pressure of 10 3 atm. This places the PET 
thruster firmly in the electrothermal 
thruster category, and justifies ignoring 
magnetic effects on the plasma acceleration 
process . 

A second effect which may be ignored is 
the compression of the plastic wall materi- 
al by the discharge pressure. For example, 
the volume compressibility of high density 
polyethylene at 10 J atra is .98, so that 
this factor reduces the electrical resis- 
tance by only 6*. Since the PFN conse- 
quently drives the current to a slightly 
higher value, the power input I Z R is nearly 
constant . 


Nozzle Processes 


For the nozzle expansion of polyethyl- 
ene (CH 2 ) plasma from To - 2 eV to T-0.5 
eV, the effective ratio of specific heats 
is about Y - 1.2, so the exit Mach number 
is 5.5 and the area r ^tio is 220. The 
pressure drops from 10^ atm to .2 atm. 
From the Saha equation we can estimate 
conditions at the exit plane. At low temp- 
eratures the degree of ionization is: 


a = 3.11 x 10 p 


,, 2 t^V‘ / 2t ,« 


For an ionization potential e - 11.2 eV for 
carbon, a =10 , so that, the average loss 
of ionization energy is .01 eV/ion, and can 
be neglected. 

Dissociation is more significant. The 
dissociation energy of the C-H bond is e - 
4.5 eV [31» and the Saha equation predicts 
a dissociation fraction p s: m j 9 
corresponding to a loss of energy in the 
exhaust stream of 2.0 eV/ion. Adding to 
the dissociation loss a thermal component 
of 5/2 T gives an approximate enthalpy in 
the exit plane of 3.3 eV/ion. For the 
assumed initial conditions the stagnation 
enthalpy is 12.6 eV/ion, so the ideal 
thruster efficiency is: 


Without an equilibrium flow nozzle, the 
plasma exhaust would be highly ionized, 
resulting in high frozen flow losses. As 
shown below in the experimental results, 
however, the ionization energy is recovered 
by three-body recombination in the nozzle. 
This rate scales aj [2]: 


(n T ) - 1 

ideal 


3.3 


T2T? 


.74 


The specific impulse for this example is 
given by I - [2 ( H -H ) ] 72 /g - 2000 
seconds. sp 0 e 


where n and T g are the electron density 
and temperature. We compare the collision- 
al recombination time 1/v~, to the particle 
flow time through the nozzle, using condi- 
tions at the exit plane where the smallest 
value of v- occurs. We find that the con- 
dition v l?<u >>>1 is easily satisfied [1], 
so that °rec-ombina t ion collisions are fre- 
quent. The flow expansion process there- 
fore remains close to Saha equilibrium, and 
because of the high expansion ratio the 
temperature drops to -0.5 eV and most of 
the ionization energy is recovered. 


Thrust Efficiency 

The high thrust efficiency of the PET 
thruster is primarily attributed to two 
factors: efficient transfer of electrical 
energy to plasma enthalpy, and efficient 
conversion or plasma enthalpy to streaming 
velocity in the nozzle. The first factor 
is a result of the high radiation absorp- 
tion coefficient or the black high density 
polyethylene propellant, and the high view 
factor or the capillary walls for the dis- 
charge black body radiation. The second 
factor is due to the ability of the high 
area ratio, high Reynolds number nozzle to 
achieve nearly equilibrium, adiabatic flow. 


For comparison we can estimate the 
performance if the flow is completely fro- 
zen at stagnation conditions. At 2 eV and 
10* 5 atm the polyethylene is -30* ionized. 
If the exhaust flow is frozen, the uissoci- 
ation and ionization energy per particle 
increases from 3.3 eV to 6.7 eV per parti- 
cle, and the efficiency drops to n - .47. 
The frozen flow specific impulse T is then 
1600 seconds. One way to view this result 
is to state that a PET thruster with an 
equilibrium flow nozzle can achieve a given 
specific impulse at a lower stagnation 
temperature and higher efficiency than a 
frozen flow device. 


Experimental Apparatus 

The relatively high impulse bit of the 
PET thruster in these experiments, typi- 
cally 0.1 N-sec, greatly simplifies the 
task or measuring thrust. The thruster, 
with mass M - 2.5 kg, is attached to a 
simple linear thrust stand (Figs. 3 and 4), 
and the recoil velocity V of several cm/sec 
is measured after each shot. The impulse 
bit is then simply /Tdt - MV[N-sec]. The 
thruster is mounted on a linear bearing and 
hardened steel shaft, and current is fed 
through flexible strips of copper mesh. A 
shorting bar in the capillary is used to 
verify zero thrust contributions from mag- 
netic fields. The recoil velocity is mea- 
sured by a linear inductive RF transducer 


[Bently Nevada Model 21504, 8 V/mm] and 
steel target, so that micron-level motions 
can be detected. The accuracy of this 
system is better than 5$. The transducer 
output signal is recorded on a digital 
oscilloscope at 2 pcec/point and is graph- 
ically differentiated to give the initial 
recoil velocity V (Fig. 5). Alternatively, 
the recoil velocity oan be determined by 
fitting the_ output signal to a waveform 
x(t) - x ° sinwt, where a, x and u can 
be easily measured (Fig. 5). °Both tech- 
niques give equivalent results. 

The thruster is sealed to the vacuum 
tank by a thin (0.4 mm) flexible diaphragm 
of rubber fabric, located at *he nozzle 
throat (Fig. 4). This arrangement keeps 
the electrical connections outside the 
vacuum tank, greatly simplifying the design 
of the high voltage insulation. The load 
on the thruster and diaphragm caused by 
atmospheric pressure is counterbalanced by 
coll springs. The natural period of the 
thruster mass-spring system is 0.090 sec, 
and the damping of the system is low (Fig. 

5) . 

Since the thrust measurement is a vital 
performance parameter a means of calibrat- 
ing the thrust stand has been devised (Fig. 

6) . A physical pendulum of mass m is con- 
structed to rotate about point 0, and im- 
pact the thruster mass m at a radius r. 
The impact point at r is designed to be the 
center of percussion about 0 so that no 
Impact forces are trmsferred to the sus- 
pension structure. The pendulum and thrust 
stand are then considered a purely rota- 
tional systenf, with angular momentum about 
0 being conserved. 

The pendulum velocity immediately be- 
fore and after Impact is found by passing 
the focussed beam of a helium-neon laser 
through a grid of equi-spaced slits (50 
line/inch Ronchi Ruling), and measuring the 
chopping frequency of the transmitted beam 
(Fig. 6). All other fixed characteristics 
of the pendulum such as the moment of 
inertia are known, so that the net momentum 
transferred to the thrust stand can be 
calculated. Simultaneously, the recoil 
momentum of the thrust stand is found using 
the position transducer (Fig. 5). Compari- 
son of the two measurements yields a cali- 
bration constant that is applied to all 
subsequent thrust data obtained using the 
transducer method. This constant is 0.985, 
indicating that the thrust stand measure- 
ment system is Inherently accurate to 
within 1.5*. 

The energy storage pulse forming net- 
work of these experiments is composed of L- 
C sections, each section consisting of a 
1.1 mH Inductor and two 29 pF high voltage 
capaoitors (Aerovox PX80D18), with a char- 
acteristic Impedance of .140 ohms. Exper- 
iments are performed with 1, 3 and 6 sec- 
tions, giving current pulse lengths of 15, 
48 and 100 pseo (Fig. 2). In all cases the 
load impedance is well matched to the 
source, resulting in efficient energy 
transfer to the discharge. 


PET thruster performance is based on 
the energy delivered to the discharge in 
the first half cycle, typically 92* of the 
stored energy, and is determined by inte- 
gration of the currrenc and voltage at the 
thruster terminals. Thruster current I is 
measured with a Rogowski loop, passively 
integrated with an RC - 2.8 msec inte- 
grator. Terminal voltage V is measured 
with a 935:1 resistive divider. Both I and 
V are recorded at 2 psei/point, enabling 
the product IV to be internally computed, 
integrated, and displayed on the scopeface 
(Fig. 7). 

The propellant is commercial high den- 
sity black polyethylene. The ablated mass 
Am is determined by weighing on a sensitive 
scale ( Voland Model 220R) to an accuracy of 
0.5 mg, and weighing again after firing IQ- 
20 shots at fixed conditions. The cathode 
and nozzle are weighed similarly. The 
nozzle ablation is quite small, and visual 
inspection after 100 shots on the nozzle 
has shown no visible erosion. The cathode 
loss is less than 15* of the total, and is 
estimated not to participate in the dis- 
charge. The capillary loss is typically 
accurate to 5*, whereas the cathode and 
nozzle erocio^ are less precise due to 
the smaller mass losses. 

Vhe thruster exhausts into a vacuum 
tank of 12 cm ID by 125 cm long, with 
access ports located every 25 cm downstream 
(Figs. 8 and 9). In addition, an access 
port is located on the centerline at the 
far end with provision for a 200 cm long 
adjustable probe mount. Several probes can 
be inserted into the tank on one shot for 
exhaust measurements, and up to three 
probes can be located at one axial down- 
stream position. 

Three tirae-ci -f 1 ight techniques have 
been employed to measure exhaust velocity: 
ion saturation probes; piezo-electric stag- 
nation pressure probes; and photodiode 
sensors (Fig. 9). The velocity measured by 
each of these techniques is generally in 
agreement. A description of each technique 
is presented below. 

The ion saturation probe is a double 
Langmuir configuration consisting of two 
identical wires, 0.23 mm diameter, embedded 
in epoxy. The wires protrude from the 
surface approximately 0.8 mm, and are spac- 
ed 1.5 mm apart. Exhaust plasma stagnates 
on the probe surface, providing a conduc- 
tive path across the wires, which are 
biased with a 29 V battery. 

In this configuration, the Langmuir 
probe signals show large, erratic fluctua- 
tions in amplitude. Suoh behavior is at- 
tributed to local cooling of the plasma due 
to evaporation of the epoxy material in the 
vicinity of the wires. However, the mean 
amplitude of the signal is consistent with 
a stagnation plasma at 2.4 eV with a den- 
sity of -io'7 oaf' 5 . Arrival timer as 
determined by the onset of the satur ion 
current signal, are in complete agreement 
with those obtained using the pressure and 
photodiode probes. 
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The stagnation pressure probe (Fig. 9) 
consists of a sensitive piezoelectric 
transducer (PCB 102A1?, 44 pai/volt), 

mounted on the vacuum tank centerline on a 
moveable support. The transducer tip 
covered by mylar tape for protection from 
the discharge plasma. Typical probe sig- 
nals are shown in Fig. 10 for two probe 
positions, showing that a t ime-of-arr i val 
can be clearly measured for each shot. 
These times are jotted and the mean slope 
is measured graphically to determine the 
exhaust velocity. 

A third t ime-of -arr i val tochnlque em- 
ploys simple red-sensitive diodes to re^'-'i 
the plasma flash on the tip of a ♦vans- 
parent plastic stinger (Fig. 9). Th probe 
is located downstream from the nozr ie exit 
plane and records a spike of radiat on when 
fast moving plasma strikes the probe (Fig. 
11). Again, the arrival times for differ- 
ent axial probe positions are plotted to 
determine the exhaust velocity. 

The total mass of the ambient gas in 
the tank at P . - .003 torr is a small 

fraction of the °exhaust mass A:a, so that 
the pressure transducer and the photodiode 
probes respond to the exhaust mass without 
interference of the background gas. The 
data from all three probes predict veloc- 
ities consistent with each other to within 
101 . 


Experimental Results 

Several series of 10 shot tests were 
run using polyethylene propellant. T he 
delivered energy ranged from 1.78 to 1.89 
kj for 15 and 48 usee pulses, a^d was 3.98 
kJ for 100 nsec pulses. The capillary 
diameter was 4.3 mm in all cases. 

The capillary length and pulse length 
were varied to permit the thruster to run 
in an unsteady or quasi-steady mode, de- 
pending on the number of acoustic travel 
times du: ing the current pulse, given oy 
ct /2l. The values of this parameter are 
gi?en in Table I below. 


Table 

1. Experimental 
Polyeth: lene 

Pulse 

Lengths for 

Test 

Pulse 

Pulse Capillary 

ct 

Series 

Energy 

, length p 

length , 

_E 


Kj 

psec 

cm 

2i 

1 . 

1.80 

15 

15 

0.4 

Unsteady 

2. 

oo 

c— 

48 

5.5 

3.5 

Quasi-steady 

3. 

1.89 

48 


4.3 

Quasi-steady 

4. 

3.98 

100 

4.5 

8.9 

Quas i -steady 


A summary of the per f oiiriAnce results is 
shown below in Table II. The estimated 
precision of the Table II measurements is: 

1. Recoil motion x(t). Impulse bit is 
measured to 5%. 

2. Exhaust probe time-of-f light , accurate 
to 10$ 

3. Current, accurate to 3$. 

4. Voltage, accurate t^ 3$. 

5. Mass loss Am, accurate to 6$. 


Di sousslon of Results 

The probe measurements of exhaust 
velocity give results which are consistent 
with our understanding of the physics of 
the discharge. In the unsteady mode cas' 
(Series 1), the current impulsively heats 
the discharge plasma which is then acceler- 
ated as the temperature decays. The mea- 
sured exhaust velocities are 20-21 km/sec. 
The quasi-steady cases (Series 2-4) achieve 
nearly steady state conditions, and hence a 
higher mean sound speed for the duration of 
the pu.se. The exhaust velocities are 
correspond ingly higher: 26- 28 km/sec. 

What evidence do we have that the ex- 
haust velocities measured correspond to 
electrothermally heated and accelerated 
propellant mass, and not to some relatively 


Table II. PET Thruster performance with Polyethylene, 10 shot averages. 


Delivered energy, E- 

1 .80 

kj 

1.78 

Maximum current , I M 

23 

kA 

17 

Pulse width, half Max 

15 

usee 

48 

Impulse bit, /Tdt 
Exhaust velocity; 

.20 

N-aec 

.095 

Pressure transducer 

21 

km/sec 

28 

Photodiode probes 

20 

km/sec 

26 

Capillary ablated mass 
Electrode ♦ nozzle 

20 

mg/shot 

6.5 

aulated mass 

NA 


NA 


kj 

i .89 

kJ 

3.98 

kj 

kA 

15 

kA 

15 

kA 

usee 

48 

usee 

100 

psec 

N-sec 

.086 

N-soc 

.189 

N-sec 

km/seo 

- 

- 

- 

- 

km/sec 

26 

km/sec 

27 

km/sec 

mg/shot 

5.3 

g/shot 

10.8 

mg/shot 


-0.9 

mg/shot 

-2.0 

mg/shot 
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low-mass precursor or shock wave? First is 
the time response of the pressure and 
photodiode probda (Figs. 10 and 11), which 
display a pulse width about equal to that 
of the current pulse. This suggests that 
the probe is responding to the momentum 
flux generated in the nozzle. Second is 
the nearly comparable velocity measured by 
the pressure and photodiode probes. Third, 
we have theoretical predictions of both the 
exhaujt velocity and probe stagnation pres- 
sure, described below, which give reason- 
ably good agreement with the measured 
values . 


Computational Model 


The predicted performance of a PET 
thruster capillary with polyethylene has 
been modeled numerically [1]. Inputs to 
the GT-Devices code include polyethylene 
chemistry [4], plasma resistivity, capil- 
lary geometry, and the experimentally mea- 
sured current and voltage waveforms. 
Because the plasma resistance is known, 
this can be used to calibrate the resis- 
tivity model used in the code, with the 
result that the plasma resistivity is close 
to that given by the Spitzer model [5]. 

The most sensitive quantity calculated 
by the modgl is th, ablated mass, which 
varies as T^ (blackbody) radiation in the 
capillary. A calculation of the ablated 
mass time history is given by the computa- 
tional model for a 2 kj shot with polyeth- 
ylene in the quasi-steady mode (Fig. 12). 
The model predicts a total ablated mass of 
6.2 mg, in close agreement with the average 
measured value of 6.5 mg per pulse. At the 
end of tie current pulse, the model pre- 
dict* that only 4.1 mg has been ejected. A 
significant fraction is apparently thus 
exhausted at a much reduced velocity, char- 
acteristic of the lower sound speed in the 
afterglow . 


2 atmospheres. The thruster stagnation 
pressure can be estimated from the measured 
impulse bit, throat area, and pulse length 
C6]: 


p ~ /Tdt ~710 atm 

V n p 

for C_ ~ 2 and t ^ 50 ysec. This pressure 
is ingood agreement with that predicted by 
the computational model (Fig. 13). If we 
assume that a normal shock sits on the 
pressure probe tip, then the probe response 
P to the shocked gas is given by [7]: 


Y 



The probe peak amplitude of 2 atm corre- 
sponds to P /P - 350. For Y - 1.2, this 
pressure ratio ^indicates M -^6.3, corre- 
sponding to a flow velocity of 25 km/sec 
for c » 9 km/sec. This value is consist- 
ent with the exhaust probe time-of-f light 
measurements . 

Based on the numerical model, the mea- 
sured exhaust velocities are consistent 
with a highly supersonic, adiabatic nozzle. 
Consistent with this model, the pulse width 
detected by the pressure transducer remains 
roughly constant with distance from the 
nozzle. This non-dispersl ve feature is 
characteristic of a high velocity puff of 
gas with a low sound speed. Furthermore, 
the measured velocity is consistent with 
adiabatic expansion of a polyethylene 
plasma with a stagnation temperature of 2.5 
eV and a stagnation pressure of -*~700 atmo- 
spheres. The stagnation enthalpy is 340 
kJ/gram for polyethylene plasma [4], which 
would then fully expand to 26 km/sec, a 
value close to that measured by the veloc- 
ity probes. 


Predicted performance curves for a 
polyethylene 2 kJ shot in the quasi-steady 
mode are shown in Fig. 13 for the capillary 
pressure, temperature, and sound speed. 
The stagnation sound speed c in the capil- 
lary is -9 krn/sec during the°current pulse, 
durirg which time Y - 1.2. Assuming adia- 
batic flow, with constant Y - 1.2 in the 
nozzle, the exit Mach number is then M - 
4.7, corresponding to the nozzle area ratio 
of 95* At this Mach number the temperature 
drops by a factor [6] 




3.2 


Performance Measurements 


The thrust to power ratio [N/kW] can be 
evaluated in the form 


T ru /,ui /Tdt /Tdt 

P [N/kW] - TTdtTfoOO Ejjtkj] 

The values of T/P for the 2 kJ shots (Table 
II) are shewn below in Table III. The 
velocities snown are those measured by the 
velocity probes for the mass exhausted 
during the current pulse. 

The thruster efficiency n T is conven- 
tionally defined as 


The sound speed thus drops by a factor of 
about /3.2 ■ 1.8, to c * 4.8 km/s. The 
corresponding predicted exhaust velocity is 
u - Me 22 km/sec, approximately the 
measured exhaust velocity. 

More evidence for high Mach number flow 
lies in the pressure probe response (Fig. 
10), which shows a maximum signal of about 


exhaust kinetic energy ^ 1 /2 mu e 
n T " energy delivered " 

For a pulsed device, the exhaust kinetic 
energy Is 

1 / 2 muf ■ / 1 / 2 illuf d t 
“ e 
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table III, Measured Thrust to Power Ratio 
for Polyethylene at 2 kj 




u ft ( Probe) 

T/P 

1 . 

Unsteady , 

15 psec pulse 

21 km/sec 

.10 

N/kW 

2. 

Quasi-steady 
48 psec pulse 

27 km/sec 

.053 

N/KW 

3. 

Quasi-steady 
48 psec pulse 

26 km/sec 

.046 

N/kW 

A. 

Quasi-steady 
100 Msec pulse 

27 km/sec 

.048 

N/kW 


where the Integral is continued past the 
end of the current pulse to account for the 
discharge afterglow, during which propel- 
lant evaporation decays to zero. Unfortu- 
nately the time dependence of m and u 
are not known, so that we oust define the 
efficiency in terns of inpulse bit as: 

_ [/Tdt] 2 

T ' 2A - E D 

where Am is the capillary ablated mass, and 
is used to define the specific impulse by 


Conclusions 

A PET thruster has been built and test- 
ed on a thrust stand, giving repeatable 
results over a series of shots. Exhaust 
velocities for polyetnylene propellant vary 
from 20 to 27 km/sec, as measured by ex- 
haust probes. These velocities are con- 
sistent with computational predictions of 
capillary conditions with an adiabatic 
expansion to high Mach number in a conical 
nozzle. Single pulse specific impulse and 
efficiency measurements based on ablated 
mass show a thruster efficiency of 37-56$ 
in the 1000-1750 second range. We believe 
that an understanding of loss mechanisms in 
the discharge would lead to an improved 
design with a thruster efficiency of 70- 
80 $. 
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Fig. 1. PET Thruster System Schematic 
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Fig. 2. PET Thruster Current Waveforms 



Fig. PET Thrust Stand Schematic 
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Fig. 5. Thrust Stand Recoil and Calibr- 
ation Response 



Fig. 3 PET Thrust Stand 



Fig. 6a. Thrust Stand Calibrator 
Schematic 
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Pig. 6b. Calibrator Optical Schematic 
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Fig. 10. Exhaust Pressure Probe Response 
at Two Positions 



Fig. 6- Vacuum Tank with Access Ports 
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